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METHOD OF CALCULATING THE FAST FOURIER TRANSFORM 
AND THE INVERSE FAST FOURIER TRANSFORM 

Abstract of the Disclosure 

A method of calculating the fast Fourier transform or the inverse fast Fourier transform 

of a series of N real samples x(n), with N power of two, operating according to a time 

interleaving algorithm and providing the sample series X(n) in ascending order to index n and 

using limited calculating and storage means. A method of calculating the fast Fourier 

transform or the inverse fast Fourier transform of a series of N conjugated complex samples 

X(n), with N power of two, operating according to a frequency interleaving algorithm. 



09/581272 

416Rec'dPCT/PTT0 0 7 JUN 2000- 

AMENDED CLAIMS 

1. A method of calculating the fast Fourier 
transform or the inverse fast Fourier transform of a 
digital signal defined by a series of N real starting 
5 samples x(n), with N a power of two and n e [0..N-1], 
comprising successive transformation steps (2) for 
transforming input samples into output samples, all the 
transformation steps being performed by means of a 
single set of butterflies with several inputs and 

10 several outputs, the operating mode of which is 
modified selectively in each transformation step, the 
input and output samples of each transformation step 
being stored in a storage memory, a series of N output 
samples y(n) representative of the fast Fourier 

15 transform or the inverse fast Fourier transform of the 
output samples x (n) being provided in the last 
transformation step, 

characterized in that output samples y(n) are 

real, 

20 and in that the output samples of a butterfly 

replace the corresponding input samples of the same 
rank in the storage memory, so that, if the starting 
samples x(n) processed in the first transformation step 
are classified in bit-reversed order of their index n, 

25 output samples y(n) are provided in the last 
transformation step in ascending order of index n, 
these output samples being defined by the following 
relations : 

y(0) = Re[X(0) ] 
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y(n) = Re[X( (n+1) /2) ] for n being odd and 
different from N-l 

y(n) = Im[X(n/2)] for- n being even and 

different from 0 

y(N-l) = Re[X(N/2) ] 
where samples X (n) , with n e [0..N-1], designate the 
complex samples of the series corresponding to the fast 
or inverse fast Fourier transform of the starting 
sample series x(n). 

2. A method of calculating the fast Fourier 
transform or the inverse fast Fourier transform of a 
digital signal defined by a series of N complex samples 
X(n) conjugated by pairs, characterized in that the 
calculation is done on a series of N real starting 
samples y(n) representative of the series of complex 
samples X(n), with N power of two and n e [0..N-1], 
the starting samples y(n) being defined as follows: 

y(0) = Re[X(0) ] 

y(n) = Re [X ( (n+1) /2) ] for n being odd and 
different from N-l 

y(n) = Im[X(n/2)] for n being even and 

different from 0 

y(N-l) = Re [X (N/2) ] 
in that this method comprises successive transformation 
steps for transforming input samples into output 
samples, a series of N real output samples x (n) 
representative of this fast or inverse fast Fourier 
transform being provided in the last transformation 
step, all the transformation steps being performed by 
means of a single set of butterflies with several 
inputs and several outputs, the operating mode of which 
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is modified selectively in each transformation step, 
the input and output samples of each transformation 
step being stored in a storage memory, 

and in that the output samples of a butterfly 
5 replace the corresponding input samples of the same 
rank in the storage memory, so that, if the starting 
samples y(n) processed in the first transformation step 
are classified in ascending order of index n, the 
output samples x(n) are provided in the last 
10 transformation step in bit-reversed order of index n. 

3. The calculation method according to claim 1 or 
2, characterized in that, in each transformation step, 
each butterfly transforms input sample pairs, the ranks 
of the input samples of the same pair within the series 

15 of input samples of said transformation step being 
symmetrical with respect to a center between the end 
rank values of the input samples transformed by said 
butterfly. 

4. The calculation method according to claim 3, 
20 characterized in that it comprises u-1 transformation 

steps E p with u=log 2 (N) and p e [0..U-2]. 

5. The calculation method according to claim 4, in 
turn dependent on claim 3, in turn dependent on claim 

25 1, characterized in further comprising: 

a preliminary step of modifying the sequence of 

the starting samples x(n) ranked in ascending order of 

index n and showing them in bit-reversed order of index 

n in the first transformation step, and 
30 - a final step of processing the series of output 

samples y(n) and providing a series of N complex 
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conjugated samples X (n) corresponding to the fast or 
the inverse fast Fourier transform of the series of 
starting samples x (n) . 

6. The calculation method of claim 4, in turn 
5 dependent on claim 3, in turn dependent on claim 1, or 
according to claim 5, characterized in that, in each 
transformation step E p , butterflies are distributed 
among N/2 P+2 calculation blocks, 

in that each calculation block has a peripheral 

10 butterfly and/or 2 P -1 internal butterflies, 

in that the peripheral butterfly of the rank a 
calculation block in transformation step Ep transforms 
the input samples of rank 2 p+2 a, 2 p+2 a+2 p+1 - 1 , 2 p+2 a+2 p+1 , 
2 p+2 a+2 p+2 -l into output samples of the same rank, 

15 and in that the internal rank x butterfly of the 

rank a calculation block in transformation step Ep 
transforms the input samples of rank 2 p+2 a+2x+l, 
2 p+2 a+2x+2, 2 p+2 a+2 p+2 -2x-3, 2 p+2 a+2 p+2 -2x-2 into output 
samples of the same rank, with p>l . 

20 7. The calculation method according to claim 6, 

characterized in that each butterfly is assigned a 
coefficient W s , whereon the calculation inside the 
butterfly is based, said coefficient being equal to e" 
j(2«s/N> with s e [0.. N/4-1] for a fast Fourier transform 

25 and is equal to e j(2l,s/N ' with s e [0.. N/4-1] for an 
inverse fast Fourier transform. 

8. Calculation method according to claim 7, 
characterized in that the internal rank x butterfly of 
the rank a calculation block in transformation step Ep 

30 is assigned coefficient W 5 with 5 = (x+1) . (N/2 P+2 ) . 
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9. The calculation method according to claim 8, 
characterized in that the butterflies for implementing 
the transformation steps are all of the same type and 
have 

5 - four inputs for receiving input samples and 

four outputs for providing output samples, 

four additional inputs, respectively primary 
mode, secondary mode, permutation, and coefficient 
inputs , 

10 in order to selectively apply different 

transformation operations to the input samples, each 
operation being determined by the values assigned to 
the primary mode, secondary mode, permutation signals, 
and a coefficient admitted on said corresponding 

15 additional inputs. 

10. The calculation method according to claim 9, 
characterized in that, for each butterfly, the primary 
mode signal is 0 for a peripheral butterfly and 1 for 
an internal butterfly, 

20 in that the permutation signal is 0 for the even 

rank calculation blocks, including rank 0, and 1 for 
the other ones. 

11. The calculation method according to claim 10, 
characterized in that, in transformation step Ep, each 

25 calculation block comprises one peripheral butterfly 
and 2 P -1 internal butterflies. 

12. The calculation method according to claim 11, 
characterized in that the secondary mode signal is 1 if 
the peripheral butterfly is used for the last 

30 transformation step, and otherwise 0. 
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13. The calculation method according to claim 12, 
characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W s =A+j.B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if the primary mode and secondary mode signals 
are 0: si = el + e2 

s2 = el - e2 
s3 = e4 - e3 
s4 = e3 + e4 

2) if the primary mode signal is 0 and the 
secondary mode signal is 1: 

si = el + e2 + e3 + e4 
s2 = el - e2 
s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

3) if the primary mode signal is 1 and the 
permutation signal is 0: 

si = el + A.e3 - B.e4 
s2 = e2 + B.e3 + A.e4 
s3 = el - A.e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 

4) if the primary mode signal is 1 and the 
permutation signal is 1: 

si = el - A.e3 + B.e4 

s2 = -e2 + B.e3 + A.e4 

s3 = el - A.e3 - B.e4 

s4 = e2 + B.e3 + A.e4 
14. The calculation method according to claim 10, 
characterized in that, in transformation step E p/ each 
calculation block comprises: 
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2 P -1 internal butterflies and a peripheral 
butterfly for the even values of index p as well as for 
the last transformation step if u is even, and 
2 P -1 internal butterflies, otherwise. 
5 15. The calculation method according to claim 13, 

characterized in that the secondary mode signal is 1 if 
the peripheral butterfly is used for the last 
transformation step with u being odd, and otherwise 0. 

16. The calculation method according to claim 15, 
10 characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W s =A+j.B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if primary mode, secondary mode and permutation 
15 signals are 0: 

si = el + e2 + e3 + e4 
s2 = el - e2 
s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 
20 2) if the primary mode signal is 0 and the 

secondary mode signal is 1: 
si = el + e4 
s2 = e2 
s3 = e3 
25 s4 = el - e4 

3) if the primary mode signal is 0 and the 
permutation signal is 1: 

si = (e3 + e4) - (el + e2) 
s2 = el - e2 
30 s3 = e4 - e3 

s4 = el + e2 + e3 + e4 
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4) if the primary mode signal is 1 and the 
permutation signal is 0: 

si - el + A.e3 - B.e4 
s2 = e2 + B.e3 + A.e4 
s3 = el - A.e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 

5) if the primary mode signal is 1 and the 
permutation signal is 1: 

si = el - A.e3 + B.e4 
s2 = -e2 + B.e3 + A.e4 
s3 = el + A.e3 - B.e4 
s4 = e2 + B.e3 + A.e4 

17. The calculation method according to claim 10, 
characterized in that, in transformation step E p , each 
calculation block comprises: 

2 P -1 internal butterflies and a peripheral 
butterfly for the even values of index p, and 
2 P -1 internal butterflies, otherwise. 

18. The calculation method according to claim 17, 
characterized in that the secondary mode signal is 1 if 
the peripheral butterfly is used for the first 
transformation step with y being even, and otherwise 0. 

19. The calculation method according to claim 18, 
characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W s =A+j.B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if the primary mode signal is 0 and the 
secondary mode signal is 1: 
si = el + e2 
s2 = el - e2 
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s3 = e4 - e3 
s4 = e3 + e4 

2) if primary mode, secondary mode and permutation 
signals are 0: 

si = el + e2 + e3 + e4 
s2 = el - e2 
s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

3) if the primary mode and secondary mode signals 
are 0 and the permutation signal is 1: 

si = (e3 + e4) - (el + e2) 

s2 = el - e2 

s3 = e4 - e3 

s4 = el + e2 + e3 + e4 

4) if the primary mode signal is 1 and the 
permutation signal is 0: 

si = el + A.e3 - B.e4 
s2 = e2 + B.e3 + A.e4 
s3 = el - A.e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 

5) if the primary mode signal is 1 and the 
permutation signal is 1: 

si = el - A.e3 + B.e4 

s2 = -e2 + B.e3 + A.e4 

s3 = el + A.e3 - B.e4 

s4 = e2 + B.e3 + A.e4 
20. The calculation method according to claim 8, 
characterized in that the butterflies for implementing 
the transformation steps are all of the same type and 
have 
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- four inputs for receiving input samples and 
four outputs for providing output samples, 

four additional inputs, respectively primary 
mode, secondary mode, permutation, and coefficient 
5 inputs, 

in order to selectively apply different 
transformation operations to the input samples, each 
operation being determined by the values assigned to 
the primary mode, secondary mode, permutation signals, 

10 and a coefficient admitted on said corresponding 
additional inputs, 

and in that the final step furthermore performs an 
addition and subtraction between the first and the last 
output sample provided in the last transformation step. 

15 21. The calculation method according to claim 20, 

characterized in that, in transformation step E p , each 
calculation block comprises one peripheral butterfly 
and 2 P -1 internal butterflies. 

22. The calculation method according to claim 21, 

20 characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W s =A+j .B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if the primary mode signal is 0: 
25 si = el + e2 

s2 = el - e2 
s3 - e4 - e3 
s4 = e3 + e4 

2) if the primary mode signal is 1 and the 
30 permutation signal is 0: 
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s2 = e2 + B.e3 + A.e4 
s3 = el - A.e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 
3) if the primary signal is 1 and the permutation 
5 signal is 1 : 

si = el - A.e3 + B.e4 
s2 = -e2 + B.e3 + A.e4 
s3 = el + A.e3 -B.e4 
s4 = e2 + B.e3 + A.e4 
10 23. The calculation method according to claim 9 or 

20, characterized in that the first and second binary- 
addresses of u bits are generated for each butterfly, 
each binary address corresponding to the rank of an 
input sample of said butterfly and the second binary 
15 address being greater than the first binary address. 

24. The calculation method according to claim 23, 
characterized in that said first and second binary 
addresses are consecutive and an internal butterfly is 
involved . 

20 25. The calculation method according to claim 23 

or 24, characterized in that, if a peripheral butterfly 
is involved, the p+2 low-order bits of the first 
address are equal to 0, and the p+2 low-order bits of 
the second address form a number equal to 2 P+1 -1 . 

25 26. The calculation method according to claim 24 

or 25, characterized in that the address of the two 
other samples to be applied to the inputs of the 
butterfly, be they peripheral or internal, are obtained 
by inverting the (p+2) low-order bits of said first and 

30 second produced addresses. 
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27. The calculation method according to claim 26, 
characterized in that even-numbered address samples and 
odd-numbered address samples are stored in two separate 
memories . 

5 28. The calculation method according to claim 25, 

characterized in that the value of the parameter s of 
the coefficient W s assigned to an internal butterfly in 
transformation step E p is coded by u-2 bits, and is: 

if p+l=p-2, the number formed by the p+1 
10 low-order bits of the second binary address produced 
for said internal butterfly, 

if p+Ku-2, the number formed by the p+1 
low-order bits of the second binary address produced 
for said internal butterfly, followed by u-p-3 zero 
15 bits at the end of the number, 

if p+l>u-2, the number formed by the p+1 low- 
order bits of the second binary address produced for 
said internal butterfly, minus its u-p-1 low-order 
bits. 

20 29. The calculation method according to claim 4, 

in turn dependent on claim 3, in turn dependent on 
claim 2, characterized in that in each transformation 
step E p/ the butterflies are distributed among 2 P 
calculation blocks, 

25 in that each calculation block comprises one 

peripheral butterfly and N/2 p+2 -l internal butterflies, 

in that the peripheral butterfly of the rank a 
calculation block in transformation step E p transforms 
the input samples of rank 2 p_p a, 2 M - p a+2 tJ - p_1 -l , 2 tJ_p a+2 p - p " 

30 l , 2 M_p a+2 iJ " p -l into output samples of the same rank, 
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and in that the internal rank x butterfly of the 
rank a calculation block in transformation step Ep 
transforms the input samples of rank 2 M ~ p a+2x+l, 2 p ~ 
p oc+2x+2, 2 p_p a+2 lJ - p " 1 -2x-3 / 2 VJ - p a+2 M_p -2x-2 into output 
samples of the same rank. 

30. The calculation method according to claim 29, 
characterized in further comprising a final step of 
modifying the sequence of the output samples provided 
in the last transformation step and classifying them in 
ascending order of index n. 

31. The calculation method according to claim 29 
or 30, characterized in that each butterfly is assigned 
a coefficient W s , whereon the calculation inside the 
butterfly is based, said coefficient being equal to e" 
j<2K S /N) with s G [0.. N/4-1] for a fast Fourier transform 
and is equal to e j(2 " s/N) with s e [0.. N/4-1] for an 
inverse fast Fourier transform. 

32. Calculation method according to claim 31, 
characterized in that the internal rank x butterfly of 
the rank a calculation block in transformation step Ep 
is assigned coefficient W 5 with 5 = (x+1) .2 P . 

33. The calculation method according to claim 32, 
characterized in that the butterflies for implementing 
the transformation steps are all of the same type and 
have 

four inputs for receiving input samples and 
four outputs for providing output samples, 

four additional inputs, respectively primary 
mode, secondary mode, permutation, and coefficient 
inputs , 
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in order to selectively apply different 
transformation operations to the input ^ samples, each 
operation being determined by the values assigned to 
the primary mode, secondary mode, permutation signals, 
5 and a coefficient admitted on said corresponding 
additional inputs. 

34. The calculation method according to claim 33, 
characterized in that, for each butterfly, the primary 
mode signal is 0 for a peripheral butterfly and 1 for 

10 an internal butterfly, 

in that the permutation signal is 0 for the even 
rank calculation blocks, including rank 0, and 1 for 
the odd values. 

35. The calculation method according to claim 31 
15 or 34, characterized in that the secondary mode signal 

is 1 if the butterfly, be it peripheral or internal, is 
used for the first transformation step, and otherwise 
0. 

36. The calculation method according to claim 35, 
20 characterized in that, for four input samples el, e2, 

e3, and e4, and for a complex coefficient W s =A+j .B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if the primary mode and secondary mode signals 
2 5 are 0: 

si = (el + e2) /2 
s2 = (el - e2) /2 
s3 = (e4 - e3) /2 
s4 = (e3 + e4)/2 
30 2) if the primary mode signal is 0 and the 

secondary mode signal is 1: 
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si = [ (el+e4) /2-e2] /2 
s2 = [ (el+e4) /2-e2] /2 
s3 - [e3- (el-e4) /2] /2 
s4 = [e3+ (el + e4) /2] /2 

3) if the primary mode signal is 1 and the 
permutation signal is 0: 

si = (el+e3)/2 
s2 = (e2+e4)/2 

S 3 = [(el-e3).A- (e2+e4).B]/2 
s4 = [-{el-e3).B + (e2+e4).A]/2 

4) if the primary mode signal is 1 and the 
permutation signal is 1: 

si = t(el-e3).A- (e2+e4).B]/2 
s2 = [-(el-e3).B + (e2+e4).A]/2 
s3 = (el+e3) /2 
s4 = (e2-e4)/2 

37. The calculation method according to claim 33, 
characterized in that the first and second binary 
addresses of u bits are generated for each butterfly, 
each binary address corresponding to the rank of an 
input sample of said butterfly and the second binary 
address being greater than the first binary address. 

38. The calculation method according to claim 37, 
characterized in that said first and second binary 
addresses are consecutive and an internal butterfly is 
involved. 

39. The calculation method according to claim 37 
or 38, characterized in that, if a peripheral butterfly 
is involved, the u-p low-order bits of the first 
address are equal to 0, and the u-p low-order bits of 
the second address form a number equal to N/2 P+1 -1. 
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40. The calculation method according to claim 38 
or 39, characterized in that the address of the two 
other samples to be applied to the inputs of the 
butterfly are obtained by inverting the u-p low-order 
bits of both produced addresses. 

41. The calculation method according to claim 40, 
characterized in that even-numbered address samples and 
odd-numbered address samples are stored in two separate 
memories. 

42. The calculation method according to claim 41, 
characterized in that the value of the parameter s of 
the coefficient W s assigned to an internal butterfly in 
transformation step E p is coded by u-2 bits, and is: 

if u-p-l=u-2, the number formed by the u-p-1 
low-order bits of the second address produced for said 
internal butterfly, 

if u-p-Ku-2, the number formed by the u-p-1 
low-order bits of the second address produced for said 
internal butterfly, followed by p-1 zero bits at the 
end of the number, 

if u-p-l>u-2, the number formed by the u-p-1 
low-order bits of the second address produced for said 
internal butterfly, minus its p+1 low-order bits. 
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METHOD OF CALCULATING THE FAST FOURIER TRANSFORM AND 
THE INVERSE FAST FOURIER TRANSFORM 

This invention relates to a method of calculating 
the fast Fourier transform or the inverse fast Fourier 
transform of a series of real numbers or a series of 
conjugated complex samples. 

The Fourier transform is probably one of the most 
important tools for analyzing, designing and 
implementing signal processing algorithms, and the 
existence of efficient algorithms, such as that of the 
fast Fourier transform, has been a major factor for 
this situation. Although most Fourier transform 
algorithms are designed for transforming series of 
complex numbers, there are, however, various 
applications, such as image or acoustic signal 
processing or certain types of multicarrier modulation 
wherein the series to be transformed are real numbers. 

In general, the direct Fourier transform and the 
inverse Fourier transform respectively set up the 
following relations between two series of N complex 
numbers, x (n) and X (n) : 

N-l 

X(n) = ^ x(k)w kn with n e [0 ... N-l] and w kn = e N 

k = 0 

x (n) = — Y X(k)w" kn with n e [0 ... N-l] 
N f-'o 

In 1965, J.W. Cooley and J.W. Tukey described an 
algorithm allowing to rapidly calculate the Fourier 
transform of a series of N complex numbers in an 
article entitled "An Algorithm for the Machine 
Calculation of Complex Fourier Series", Math. 



Computation, Vol. 19, 1965, pp. 297-301. This algorithm 
is very interesting when N is a power of two because it 
is then particularly straightforward to implement. This 
algorithm requires u calculation steps where u=log 2 (N) . 
It is based on breaking down the series to be 
transformed into two interleaved sub-series. There are 
two kinds of interleaving: time interleaving and 
frequency interleaving. Both kinds of interleaving are 
explained more in detail in the course of the 
description with reference to Figs. 1A and IB. 

Specific algorithms have been developed from this 
basic algorithm to deal with the case of real number 
series. The Fourier transformation of a series of 8 
real numbers according to a time interleaving algorithm 
and according to a frequency interleaving algorithm is 
illustrated in Figs. 1A and IB. For each time 
interleaving Fourier transform algorithm there is a 
frequency interleaving algorithm, which corresponds to 
a double inversion of the series of transformation 
operations, on the one hand, and for each butterfly 
circuit, of the proper transformation operations, on 
the other hand. Whatever the interleaving chosen, the 
transformation method requires three transformation 
steps E0, El, and E2, these steps being implemented 
through a set of four butterfly circuits CC, generally 
known as "butterflies" in technical speak. Each 
butterfly, represented in the figure by a point of 
intersection between two columns of numbers, performs 
calculations on two numbers, real or complex. The 
symbols R and C respectively identify a real number and 
a complex number. The sequence of real and/or complex 



numbers after the transformation steps depends on the 
interleaving chosen. 

The time interleaving algorithm is generally- 
chosen for calculating the Fourier transform of a 
series of real numbers because of the symmetrical 
distribution of real and complex numbers throughout the 
steps. On the other hand, the frequency interleaving 
algorithm is more suitable for the direct or inverse 
transformation of a series of conjugated complex 
numbers . 

When the series to be transformed x(n) is real, 
the Fourier transform verifies the following relation: 

x (n) is real if and only if 

X(n) = X* (-n) =X* (N-n) ; (1) 
where * designates the conjugating operation. 

For a series x(n) of N real numbers, the following 
results are inferred from this relation: 

- X(0) and X(N/2) are real; 

- X(n) = X* (N-n) for 1 < n < N/2 -1 

Relation (1) highlights the presence of redundant 
information in the X(n) series. 

It should be noted that the transformation method 
is generally implemented by a single set of 
butterflies, the operating mode of which is modified as 
the transformation goes along. At each change of 
operating mode, results are stored in a memory having N 
storage locations, the output samples of a butterfly 
replacing the corresponding input samples of the same 
rank in the memory. This method of applying the 
algorithm is generally known as the "in place" method. 
This method has a major advantage: if the elements of 
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the x(rt) series are processed in the first 
transformation step in bit-reversed order of index n, 
the numbers of the X{n) series are output in the last 
transformation step in ascending order of index n and 
5 vice versa. 

A known transformation method is shown by way of 
example in Fig. 2. This method performs the Fourier 
transformation of a real x(n) series according to a 
complex time interleaving algorithm. In this example, 

10 the x (n) series to be transformed comprises sixteen 
real samples, x(0) to x(15). The transformation method 
comprises four transformation steps E p with 0<p<3 . The 
samples of the x(n) series are shown in the first 
transformation step in bit-reversed order of their 

15 index n . 

At this stage of the explanations, the terms used 
in the course of the description should be defined. The 
rank of a sample is taken from the position it occupies 
in the series of samples to which it belongs. The index 
20 of a sample then corresponds to the starting rank of 
this sample. 

The intermediate results obtained in the various 
transformation steps are represented by the series 
A(n) , B(n), and C(n). The samples of the series x(n), 

25 A(n), B(n), C(n), and X (n) are stored in double storage 
locations, one storage location being reserved for the 
real portion of the sample and the other location being 
reserved for the imaginary portion thereof. A R (n) and 
A x (n) respectively designate the real portion and the 

30 imaginary portion of the index n sample of the A(n) 
series. Butterflies are represented in the figure by 
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points of intersection between columns of storage 
locations. Each butterfly is assigned a coefficient W s 
symbolized in Fig. 2 by a pair of coordinates A/B where 
A and B respectively designate the real portion and the 
5 imaginary portion of the coefficient W s . Coordinates 
1/0 and 0/-1 are respectively assigned to coefficients 
W°=l and W N / 4 =W 4 =-j . For the sake of clarity and in 
order to simplify their formulation, the remaining 
coefficients W s have been represented by the following 
10 pairs: 

W 1 -> 2 / -4 W 5 -> -4 / -2 

W 2 -> 3 / -3 W 6 -> -3 / -3 

W 3 -> 4 / -2 W 7 -> -2 / -4 

These pairs of coordinates are graphically 

15 represented in Fig. 3. In fact, coordinates A and B 
respectively represent a cosine value and a sine value. 
This coefficient W s participates in the calculation 
performed by the butterfly. Furthermore, the 
butterflies are distributed at each transformation step 

20 among N/2 P+1 calculation blocks, each calculation block 
comprising 2 P butterflies. In the course of the 
description, the parameter q designates the rank of the 
calculation blocks within the same transformation step; 
q is included between 0 and (N/2 P+1 )-1. 

25 In the first transformation step E 0/ the 

butterflies are distributed among eight calculation 
blocks, each comprising a butterfly performing an 
operation on two complex or real samples. If el and e2 
are to designate the samples applied to the inputs of a 

30 butterfly, the latter outputs samples si and S2 defined 
as follows: 
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si = el + W s .e2 and s2 = el -(W s .e2) 
where W s is the coefficient assigned to said butterfly. 

For this first transformation step, coefficient 
W°=l is assigned to the eight butterflies. As the 
5 samples x (n) and the coefficient W° are real, the 
samples A(n) obtained at the end of step E 0 are real. 

For the second transformation step, Ei, the 
butterflies are distributed among four calculation 
blocks each comprising two butterflies. Coefficient 
10 W°=l is assigned to the first one of these butterflies; 
thus, the first butterfly of each calculation block 
provides two real samples. The second butterfly of the 
calculation blocks is associated with coefficient 
W N/4 =W 4 =-j and generates two conjugated complex samples. 
15 The output samples obtained at the end of step E x are 
designated by the series B(n). 

For the third transformation step, E 2 , the 
butterflies are distributed among two calculation 
blocks each comprising four butterflies respectively 
20 associated with coefficients W°, W 2 , W 4 , and W 6 . The 
output samples in step E 2 are designated by the series 
C(n) . Finally, for the fourth transformation step, E 3 , 
a single calculation block comprising eight butterflies 
respectively associated with coefficients W°, W 1 , W 2 , 
25 W 3 , W 4 , W 5 , W 6 , and W 7 , is provided. This transformation 
step generates the transformed series X (n) . 

Given relation (1), the X (n) series comprises, on 
the one hand, real samples, X(0) and X(8), and on the 
other hand, complex samples, X(l) to X(7) and X(9) to 
30 X(15), samples X(15) to X(9) respectively being the 
conjugates of samples X(l) to X(7). The X (n) series 



therefore contains redundant information. The storage 
locations outlined in thick stroke in Fig. 2 designate 
the storage locations enclosing the conjugated values 
of the complex samples contained in the storage 
locations associated therewith by an arrow. The 
intermediate results series B(n) and C(n) also contain 
redundant information. 

It is then possible to delete this redundant 
information in order to reduce by half the size of the 
sample storage memory as well as the number of 
butterflies . 

However, removing redundant information stored in 
the storage locations outline in bold strokes in Fig. 2 
implies the complete reorganization of the 
transformation steps of Fig. 2. Thus reorganizing the 
transformation has the effect of modifying the output 
sequence of samples X (n) . 

The problem therefore consists in reducing the 
size of the storage memory and the number of 
butterflies while maintaining the output sequence of 
the X(h) samples. It is the object of the invention to 
offer a method of calculating the fast Fourier 
transform or the inverse fast Fourier transform of a 
series of N real samples x(n), with N power of 2, 
operating according to a time interleaving algorithm, 
which provides the series of samples X (n) in ascending 
order of index n and uses limited calculation and 
storage means . 

For this purpose, the object of the invention is a 
method of calculating the fast Fourier transform or the 
inverse fast Fourier transform of a digital signal 



defined by a series of N real starting samples x (n) , 
with N power of two and n e [0..N-1], comprising 
successive transformation steps for transforming input 
samples into output samples, all the transformation 
5 steps being performed by means of a single set of 
butterflies with several inputs and several outputs, 
the operating mode of which is modified selectively in 
each transformation step, the input and output samples 
of each transformation step being stored in a storage 
10 memory, a series of N output samples y(n) 
representative of the fast Fourier transform or the 
inverse fast Fourier transform of the starting samples 
x(n) being provided in the last transformation step, 

characterized in that output samples y(n) are 

15 real, 

and in that the output samples of a butterfly 
replace the corresponding input samples of the same 
rank in the storage memory, so that, if the starting 
samples x (n) processed in the first transformation step 
20 are classified in bit-reversed order of their index n, 
output samples y(n) are provided in the last 
transformation step in ascending order of index n, 
these output samples being defined by the following 
relations : 
25 y(0) = Re[X(0) ] 

y(n) = Re [X ( (n+1 ) / 2) ] for n being odd and 

different from N-l 
y(n) = Im[X(n/2)] for n being even and 

different from 0 

30 y (N-l) = Re [X (N/2) ] 
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where the X(n) samples, with n e [0..N-1] designate 
the complex samples of the series corresponding to the 
fast or inverse fast Fourier transform of the series of 
starting samples x (n) . 
5 For methods operating according to a frequency 

interleaving algorithm, the invention also relates to a 
method of calculating the fast Fourier transform or the 
inverse fast Fourier transform of a digital signal 
defined by a series of N complex samples X (n) 
10 conjugated by pairs represented by a series of N real 
starting samples y(n), with N power of two and n e 
[0..N-1], the starting samples y(n) being defined as 
follows: 

y(0) = Re[X(0) ] 

15 y(n) - Re [X ( (n+1) /2) ] for n being odd and 

different from N-l 
y(n) = Im[X(n/2)] for n being even and 

different from 0 

y (N-l) - Re[X (N/2) ] 

20 this calculation method comprising successive 
transformation steps for transforming input samples 
into output samples, a series of N output samples x (n) 
representative of this fast or inverse fast Fourier 
transform being provided in the last transformation 

25 step, all the transformation steps being performed by 
means of a single set of butterflies with several 
inputs and several outputs, the operating mode of which 
is modified selectively in each transformation step, 
the input and output samples of each transformation 

30 step being stored in a storage memory, 
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characterized in that output samples x (n) are 

real, 

and in that the output samples of a butterfly 
replace the corresponding input samples of the same 
5 rank in the storage memory, so that, if the starting 
samples y(n) processed in the first transformation step 
are classified in ascending order of index n, the 
output samples x(n) are output in the last 
transformation step in bit-reversed order of index n. 

10 The inventive calculation methods perform 

operations on real samples and consequently use limited 
calculation and storage means in comparison with the 
method shown in Fig. 2. 

According to another feature of the invention, in 

15 each transformation step, the butterflies transform 
input sample pairs, the ranks of the input samples of 
the same pair within the series of input samples of 
said transformation step being symmetrical with respect 
to a center between the end rank values of the input 

20 samples transformed by said butterfly. Input samples 
processed by the same butterfly are thus symmetrically 
linked together by pairs. The result is simplified 
handling of sample addressing. 

According to another aspect of the invention, the 

25 method preferably comprises u-1 transformation steps E p 
with u=log 2 (N) and p e [0..U-2]. 

Other features and advantages of the invention 
will be apparent from reading the following detailed 
description, which is made with reference to the 

30 appended drawings, where: 
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- Figs. 1A and IB, already described, 
respectively represent a Fourier transformation of 
eight real numbers according to a time interleaving 
algorithm and according to a frequency interleaving 

5 algorithm; 

Fig. 2, already described, illustrates the 
transformation of a series of 16 real numbers into a 
series of 16 complex numbers according to a complex 
time interleaving algorithm; 
10 - Fig. 3, already described, graphically 

represents the mapping of coefficients W s and 
coordinate pairs A/B; 

- Fig. 4 represents modifications applied to part 
of the transformation of Fig. 2; 

15 - Fig. 5 illustrates a modified transformation 

only processing real numbers; 

Fig. 6 is a representation of the method of 
calculating the fast Fourier transform according to the 
invention; 

20 - Figs. 7A and 7B respectively illustrate 

permutations performed on peripheral butterflies and on 
internal butterflies with odd-numbered rank of the 
transformation of Fig. 5; 

Fig. 8 represents an embodiment of the 

25 transformation method according to the invention, 
comprising u transformation steps; 

Fig. 9 illustrates a grouping of peripheral 
butterflies according to a first embodiment of a 
transformation method comprising u-1 transformation 

30 steps; 
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- Fig. 10 represents a first embodiment of a 
transformation method comprising u-1 transformation 
steps ; 

Fig. 11 represents a butterfly design relating 
5 to the transformation method illustrated in Fig. 10; 

- Fig. 12 represents an alternative of the 
embodiment of Fig. 10; 

Fig. 13 represents a butterfly design relating 
to the embodiment of Fig. 12; 
10 - Fig. 14 illustrates a grouping of peripheral 

butterflies according to a second embodiment of a 
transformation method comprising u-1 transformation 
steps ; 

Fig. 15 represents a second embodiment of a 
15 transformation method comprising u-1 transformation 
steps, with p being even; 

Fig. 16 represents an alternative of the 
preceding embodiment, with u being odd; 

- Fig. 17 represents a butterfly design relating 
2 0 to the embodiments shown in Figs. 15 and 16; 

- Fig. 18 represents a third embodiment of a 
transformation method comprising u-1 transformation 
steps; 

Fig. 19 represents a butterfly design relating 
25 to the embodiment of Fig. 18; 

Fig. 20 represents the addresses that are 
associated with the various butterflies implemented in 
the embodiment shown in Fig. 12; 

- Fig. 21 represents the addresses that are 
30 associated with part of the samples of a transformation 

method processing a series of 32 real samples; 
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Fig. 22 represents a sample embodiment of a 
transformation method operating according to a 
frequency interleaving algorithm; 

Fig. 23 represents a butterfly design relating 
5 to the embodiment of Fig. 22; 

Fig. 24 represents the addresses that are 
associated with the various butterflies of the 
embodiment illustrated in Fig. 22. 

According to the invention, only one part of the 

10 X(n) samples is calculated, the other part of the 
samples being redundant. E.g., the calculation could be 
limited to the X(n) numbers for 0<n<N/2. However, a 
more adequate solution is calculating the X(n) numbers 
for 0<n<N/4 and N/2<n<3N/4. The latter solution is 

15 preferable because it only involves the first (N/4)+l 
butterflies for calculating X(n), as the last (N/4)-l 
butterflies can be removed. The storage locations thus 
released can be used for storing the real portion or 
the imaginary portion of the remaining X (n) numbers. 

20 The size of the storage memory can thus be halved if 
the size of the storage locations is limited to storing 
a real number instead of a complex number. The real 
portion of the complex X (n) number is stored in the 
storage location that has been assigned thereto 

25 initially whereas the imaginary portion is stored in 
the storage location initially assigned to the number 
X (N-n) . The same operation can be performed for the 
series of intermediate results B(n) and C(n) . The 
transformation method of the invention is limited to 

30 calculating N real samples y(n) instead of N complex 
samples. The real samples y(n) are defined as follows: 
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y(0) = Re[X(0) ] 

y(n) = Re [X ( (n+1) /2) ] for n being odd and 

different from N-l 
y(n) = Im[X(n/2) ] for n being even and 

5 different from 0 

y(N-l) - Re[X(N/2) ] . 

Removing redundant information and reorganizing 
storage locations deeply modifies the diagram of 
Fig. 2. Fig. 4 illustrates the reorganizations 

10 performed on the part relating to the calculation of 
samples A(l), A(3), A(5), A(7), B(l), B(3), B (5) , B(7), 
C(l), C(3), C(5), and C(7). Redundant samples to be 
removed are B(3), B(7), C(5), and C(7). When the size 
of the storage locations has been reduced and the 

15 samples remaining in theses storage locations have been 
reorganized, the butterflies perform calculations on 
real numbers. The butterflies to which the pair 0/-1 is 
assigned perform calculations on two real numbers. In 
practice, they copy onto their first output the number 

20 present at their first input and multiply by -1 the 
number present at their second input and provide it to 
their second output. The butterflies to which the pair 
1/0 is assigned perform an addition and a subtraction 
on two real numbers. Finally, the other butterflies 

25 perform operations on four real numbers. 

The transformation thus rearranged is illustrated 
in Fig. 5. In this figure, the butterflies associated 
with the pairs 1/0 and 0/-1 corresponding to the 
coefficients W° and W 4 are connected with storage 

30 locations by thick strokes. This figure shows that the 
reorganization of the transformation method steps 



modifies the output sequence of the y(n) samples and 
therefore the output sequence of the X (n) series. 
Furthermore, this transformation method no longer has a 
specific symmetry allowing to link by pairs the samples 
processed by the same butterfly. The result is a very 
complicated address management of the samples to be 
applied to the butterfly inputs. 

Fig. 6 shows the steps of a method of calculating 
the fast Fourier transform or the inverse fast Fourier 
transform of a series of N real numbers x(n), with N 
power of 2, operating according to a time interleaving 
algorithm. It mainly comprises transformation steps 2 
for transforming N starting samples x(n) classified in 
bit-reversed order of their index n into real output 
samples y(n) representative of this Fourier transform 
classified in ascending order of index n. 

Advantageously, it comprises a preliminary step 1 
for ranking the N real starting samples x(n) to be 
transformed in the bit-reversed order of their index n 
if the samples x (n) are not already in this sequence 
and a final step 3 for generating the N complex samples 
X(n) corresponding to the fast Fourier transform of the 
starting samples x (n) from the N real samples y(n) 
obtained at the end of the transformation steps. 

The methods that will be detailed in the course of 
the description will be more in particular for 
calculating the fast Fourier transform of a real 
series. Also, the coefficients W s assigned to the 
butterflies for implementing the inventive method will 
be of the type e -j< 2 " s/N > with s>0 . For calculating the 
inverse fast Fourier transform the calculation method 



is the same, however, the coefficient is of the type 
e j(2«s/N) with s > 0 _ 

In order to obtain at the same time y(n) samples 
sorted in ascending order of index n and symmetry of 
calculation, according to the invention, it is 
suggested to modify the calculations performed by the 
butterflies of the odd rank g calculation blocks of the 
transformation illustrated in Fig. 5 in accordance with 
the diagrams of Figs. 7A and 7B. 

As the butterflies associated with coefficient 1/0 
(Fig. 7A) of the odd rank calculation blocks are 
involved, provision is made for permutation both 
outputs of the butterfly and multiplying by -1 the 
result provided at the second butterfly output. 

As the four input butterflies (Fig. 7B) are 
involved, provision is made for permutation the first 
two outputs with the last two ones. 

This method is applied to the whole transformation 
and then a method is obtained providing at the output 
y(n) samples in ascending order of index n. This method 
is illustrated in Fig. 8. The symbol^ placed above the 
calculation blocks designates the calculation blocks 
wherein butterflies have been modified, i.e. odd rank 
calculation blocks. Due to the symmetry of calculation 
of the transformation method, the intersecting points 
representing the butterflies are superposed inside each 
calculation block. 

The butterflies to which the coefficient 1/0 is 
assigned are called peripheral butterflies because they 
perform calculations on samples arranged at the ends of 
the calculation block. The other butterflies are called 
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internal butterflies. It should be noted that in each 
transformation step, not all the samples are processed 
always and that unprocessed samples are kept in their 
storage locations to be processed in subsequent steps, 
5 or else produced as output if they are already in their 
final shape. 

The transformation method thus modified provides 
y(n) samples in ascending order of index n and, in each 
transformation step, has symmetry of calculation 
10 facilitating the addressing of the samples to be 
processed . 

According to another aspect of the invention, the 
calculation method advantageously comprises u-1 
transformation steps. Several embodiments derived from 
15 that of Fig. 8 and comprising u-1 transformation steps 
are therefore shown in the course of the description. 
One butterfly design is associated with each of these 
embodiments. 

All of these embodiments have the following 
2 0 features in common: 

in each transformation step, provision is made 
for N/2 P+2 calculation blocks and each calculation block 
comprises an peripheral butterfly and/or 2 P -1 internal 
butterflies; all butterflies, be they peripheral or 
25 internal, perform calculations on four real samples; 

- the ranks of the samples processed by the same 
butterfly are defined as follows: if, in transformation 
step Ep, a peripheral butterfly belonging to the rank a 
calculation block is taken into consideration, it 
30 transforms the input samples of rank 2 p+2 a, 2 p+2 a+ 2 P+1 - 1 , 
2 p+2 a+2 p+1 , 2 p+2 a+2 p+2 -l into output samples of the same 
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rank, and, if an internal rank x butterfly in a rank a 
calculation block is taken into consideration in step 
Ep, it transforms the input samples of rank 2 p+2 a+2x+l, 
2 p+2 a+2x+2, 2 p+2 a+2 p+2 -2x-3 / 2 p+2 a+2 p+2 -2x-2 into output 
samples of the same rank, with P>1; 

the coefficient assigned to the internal rank x 
butterfly of the rank a calculation block in step Ep is 
equal to W 5 with 5 = (x+1) . (N/2 P+2 ) . 

In all of these embodiments, the input samples of 
each butterfly form sample pairs, the sample ranks of 
the same pair in the series of input samples of a 
transformation step being symmetrical with respect to 
the center value of the end ranks of the input samples 
transformed by said butterfly. This center corresponds 
to the value 2 p+2 a+2 p+1 -l/2 . Therefore, one just has to 
know the rank of two of the four samples to be applied 
to the inputs of the butterfly in order to infer 
therefrom the rank of the other two. Addressing these 
samples is thus simplified. This will be explained more 
in detail at a later point of the description. 

Thus, according to a first embodiment, the 
neighboring calculation blocks in each transformation 
step are grouped by pairs. The peripheral butterflies 
of the same calculation block are then merged into a 
single peripheral butterfly. A sample merging of two 
peripheral butterflies is shown in Fig. 9. This example 
relates to peripheral butterflies associated with 
samples x(0), x(8), x(4), and x (12) . 

Furthermore, as the transformation steps do not 
process all samples each time, certain calculations can 



be anticipated. E.g., calculating samples C R ( 1 ) , Ci ( 1 ) , 
C R (3), and Ci(3) can be done in the second 
transformation step. The result is the diagram of Fig. 
10 showing a first embodiment of the transformation 
circuit wherein the method only comprises u-1 
transformation steps. In transformation step E p , each 
calculation block has a peripheral butterfly and 2 P -1 
internal butterflies. It can be considered that this 
method only comprises 3 transformation steps, the 
fourth step being limited to performing an addition and 
a subtraction. This addition and this subtraction are 
preferably performed during the final step, and in 
order to limit the number of transformation steps. 

A butterfly design associated with the embodiment 
of Fig. 10 is represented in Fig. 11. It comprises: 

four inputs for receiving input samples el, e2, 
e3, e4, and four outputs for providing output samples 
si, s2, s3, s4, and 

three additional, respectively primary mode MP, 
permutation PERM, and coefficient COEF, inputs. 

This butterfly is responsible for selectively 
applying to input samples el, e2, e3, and e4, various 
transformation operations each determined by the values 
assigned to primary mode, permutation signals and to 
coefficient W s admitted at the corresponding additional 
inputs . 

The primary mode signal is 0 for a peripheral 
butterfly and 1 for an internal butterfly. When the 
permutation signal is a 1, the output samples si and s2 
of the butterfly are swapped with output samples s3 and 
s4. This permutation is only possible if the butterfly 
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is an internal one. Finally, the coefficient W 
associated with the butterfly is applied to the 
coefficient input COEF. 

Thus, if the complex coefficient W s =A+j.B is 
applied to the coefficient entry of the butterfly, the 
latter provides the following output samples si, s2, 
s3, and s4 

1) if the primary mode signal is 0: 
si = el + e2 

s2 = el - e2 
s3 = e4 - e3 
s4 = e3 + e4 

2) if the primary mode signal is 1 and the 
permutation signal is 0: 

si = el + A.e3 - B.e4 
s2 = e2 + B.e3 + A.e4 
s3 = el - A.e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 

3) if the primary signal is 1 and the permutation 
signal is 1 : 

si = el - A.e3 + B.e4 

s2 = -e2 + B.e3 + A.e4 

s3 = el + A.e3 -B.e4 

s4 = e2 + B.e3 + A.e4 
According to an alternative embodiment, provision 
can be made for addition and subtraction to be 
performed inside the peripheral butterfly of the last 
transformation step. This embodiment is shown in Fig. 
12. For this purpose, the corresponding butterfly 
design has a fourth additional input called secondary 
mode input MS to which a secondary mode signal is 
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applied. This signal is 1 for the peripheral butterfly 
of the last transformation step, otherwise it is 0. 
This design is illustrated in Fig. 13. This design has 
an additional operating mode in comparison with the 
5 preceding one; thus, when the primary mode signal is 0 
and the secondary mode signal is 1, the output obtained 
is : 

si = el + e2 + e3 + e4 
s2 = el - e2 
10 s3 = e4 - e3 

s4 - (el + e2) - (e3 + e4) 

According to a second embodiment derived from the 
diagram of Fig. 8, the neighboring peripheral 
butterflies in even index p transformation steps E p are 

15 grouped by pairs and are merged with the peripheral 
butterfly of the second odd index step in order to form 
a new peripheral butterfly at the odd index step. This 
grouping is illustrated in Fig. 14 in an example. In 
this example, the peripheral butterflies of the first 

20 transformation step processing samples x(0), x(8), 
x(4), and x(12) are merged with the peripheral 
butterfly of the second step processing samples A R (0) 
and A R (2). Also, the peripheral butterflies of the 
first transformation step processing samples x(2), 

25 x(10), x(6), and x(14) are merged with the peripheral 
butterfly of the second step processing samples A R ( 4 ) 
and A R (6). The two butterflies obtained are different 
in that the second one performs in addition a 
permutation between the first and second outputs. If 

30 this grouping is applied to the whole transformation 
illustrated in Fig. 8, the result is that step E 0 no 
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longer uses any butterflies and can be removed. The 
resulting transformation method is illustrated in Fig. 
15. 

However, two cases should be distinguished for 
5 this transformation method: the case where N is an even 
power of two (u being even) and the case where N is an 
odd power of two (u being odd) . 

In case u is even, there is an even number of 
transformation steps in the embodiment shown in Fig. 8 
10 and grouping the peripheral butterflies of even index 
steps with those of the following odd index steps is no 
problem. This case corresponds to the diagram in Fig. 
15. 

In case u is odd, the peripheral butterfly of the 

15 last even index step cannot be grouped with other 
peripheral butterflies. Therefore, a specific operating 
mode should be provided for this case. This case is 
illustrated in Fig. 16, this figure representing the 
transformation of a series of eight real samples (u=3) . 

20 The peripheral butterfly of the last transformation 
step of this circuit could not be merged with other 
peripheral butterflies. 

The butterfly design associated with this second 
embodiment is illustrated in Fig. 17; it differs from 

25 the preceding design in that the secondary mode signal 
is 1 when a peripheral butterfly for implementing the 
last step is involved and u is odd, and in that 
permutation applies to all the butterflies of the even 
rank calculation blocks. 

30 The calculations performed by the butterfly are 

also different and are defined as follows: 
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1) if primary mode, secondary mode and permutation 
signals are 0: 

si = el + e2 + e3 + e4 
s2 = el - e2 
5 s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

2) if the primary mode signal is 0 and the 
secondary mode signal is 1: 

si = el + e4 
10 s2 = e2 

s3 = e3 
s4 = el - e4 

3) if the primary mode signal is 0 and the 
permutation signal is 1: 

15 si = (e3 + e4) - (el + e2) 

s2 = el - e2 
s3 - e4 - e3 
s4 = el + e2 + e3 + e4 

4) if the primary mode signal is 1 and the 
20 permutation signal is 0: 

si = el + A.e3 - B.e4 

s2 = e2 + B.e3 + A.e4 

s3 = el - A.e3 + B.e4 

s4 = -e2 + B.e3 + A.e4 
25 5) if the primary mode signal is 1 and the 

permutation signal is 1: 

si = el - A. e3 + B . e4 

s2 = -e2 + B.e3 + A.e4 

s3 = el + A.e3 - B.e4 
30 s4 = e2 + B.e3 + A.e4 
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In case u is odd, it is also possible on the one 
hand to provide for the peripheral butterflies to be 
grouped for implementing the first transformation step 
in the same way as in the first embodiment, and on the 
5 other hand, for the butterflies of the other steps to 
be grouped as in the third embodiment. Grouping 
butterflies from the second step on is then performed 
by assigning an even index to the first transformation 
step. These groupings are represented in Fig. 18. 

10 The butterfly design corresponding to this 

embodiment is represented in Fig. 19. The secondary 
mode signal is 1 for a peripheral butterfly 
implementing the first transformation step of the 
circuit and if u is even. The calculations performed by 

15 this butterfly are the following ones: 

1) if the primary mode signal is 0 and the 
secondary mode signal is 1: 

si = el + e2 
s2 = el - e2 
20 s3 = e4 - e3 

s4 = e3 + e4 

2) if primary mode, secondary mode and permutation 
signals are 0: 

si = el + e2 + e3 + e4 
25 s2 = el - e2 

s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

3) if the primary mode and secondary mode signals 
are 0 and the permutation signal is 1 : 

30 si - (e3 + e4) - (el + e2) 
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s3 = e4 - e3 

s4 = el + e2 + e3 + e4 

4) if the primary mode signal is 1 and the 
permutation signal is 0: 

5 si = el + A.e3 - B.e4 

s2 = e2 + B.e3 + A.e4 
s3 = el - A.e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 

5) if the primary mode signal is 1 and the 
10 permutation signal is 1: 

si = el - A.e3 + B.e4 
s2 = -e2 + B.e3 + A.e4 
s3 = el + A.e3 - B.e4 
s4 = e2 + B.e3 + A.e4 
15 In all the embodiments described before, the ranks 

of samples processed by the same butterfly are 
symmetrical by pairs with respect to a center value. 
One just has to know the rank of the first two input 
samples of the butterfly to infer the other two ones 
20 therefrom by symmetry. If the input and output samples 
are saved in storage locations the address of which 
corresponds to the rank of these samples, addressing 
the latter will be simplified. 

Indeed, all that needs to be done then is to 
25 generate two addresses per butterfly, as the other two 
ones can be inferred by symmetry. Furthermore,- it 
should be noted that the addresses of the input samples 
and those of the corresponding output samples are the 
same because the transformation is applied according to 
30 an in place method. 



The addresses associated with the various 
butterflies relating to the transformation method of 
Fig. 12 are shown in Fig. 20. The address of a sample 
is taken to be equal to the rank of this sample in the 
5 sample series to which it belongs. The series to be 
transformed in the example of Fig. 12 comprises 16 
samples. Therefore, 16 addresses have to be produced, 
with four-bit addressing. For a series of N samples, 
log 2 (N) -bit addressing is performed. 

10 Each transformation step of the method of Fig. 12 

is implemented by four butterflies each processing four 
real input samples. The binary addresses of the samples 
to be processed in each step are therefore distributed 
among four groups of four addresses. The address groups 

15 relating to one peripheral butterfly are contained in 
bold line boxes and the address groups relating to an 
internal butterfly are contained in thin line boxes. 
Furthermore, the address groups relating to the same 
calculation block are grouped in dotted line boxes. 

20 The first peripheral butterfly in step E 0 

processes the first four samples of the x(n) series 
resulting from the preliminary classification step. The 
corresponding addresses to be generated for this 
butterfly are therefore 0000, 0001, 0010, and 0011. 

25 Also, considering the first peripheral butterfly of 
step Ei, it processes the samples of rank 0, 3, 4, and 
7 of the A(n) series. The corresponding addresses to be 
generated for this butterfly are therefore 0000, 0011, 
0100, and 0111. 

30 The (u-2-p) high-order bits of the addresses 

indicate rank q of the calculation block to which the 
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address is related. Thus, the two high-order bits of 
the addresses generated for the first calculation block 
of step E 0 are 00. It should also be noted that, as the 
last step only uses a single calculation block, the 
5 addresses generated for this calculation block have no 
bit relating to the rank of this block (u-p-2=0) . 

For each butterfly, provision is made for 
generating only two binary addresses; the other two 
ones are obtained by inverting the (p+2) low-order bits 

10 of the generated addresses. E.g., considering the 
addresses of the first peripheral butterfly in step E x , 
only addresses 0000 and 0011 are generated and 
addresses 0111 and 0100 are obtained by inverting the 3 
low-order bits of the generated addresses. 

15 A first and a second address per butterfly are 

then produced, these addresses being consecutive for an 
internal butterfly. When a peripheral butterfly is 
involved, the p+2 low-order bits of the first address 
are equal to 0, and the p+2 low-order bits of the 

20 second address form a number equal to 2 P+1 -1. 

Regarding the coefficients W s to be applied to the 
coefficient inputs COEF of the internal butterflies, 
they are stored in a memory of the calculation circuit. 
Only N/4 coefficient values are required for 

25 calculating the Fourier transform, u-2-bit addressing 
is performed for these coefficients. In the embodiments 
described previously where N=16, only the pairs 1/0, 
2/-4, 3/-3, and 4/-2 are used which correspond to 
coefficients W°, W 1 , W 2 , W 3 . 

30 According to the invention, the address of these 

four coefficients must therefore be known. The address 
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associated with each coefficient W s is chosen to be 
equal to the value of the power s. Consequently, the 
addresses of coefficients W°, W 1 , W 2 , W 3 are 00, 01, 10, 
11, respectively. 
5 So as not to be obliged to generate these 

addresses, according to the invention, the addresses 
produced for addressing the samples are used. The 
address of the coefficient which is assigned to a 
butterfly is included in the second address produced 
10 corresponding to the greater one of both addresses. 

However, we must distinguish between three cases: 

a) when p+l=u-2, the coefficient address 
corresponds to the number formed by the p+1 low-order 
bits of the second address generated for this internal 

15 butterfly. This is the case for the second step (p=l) 
in the example of Fig. 20. The 2 low-order bits of the 
second address are 10 and therefore designate 
coefficient W 2 . 

b) if p+l>u-2, the coefficient address corresponds 
20 to the number formed by the p+1 low-order bits of the 

second address generated for this internal butterfly, 
minus its u-p-1 low-order bits. This is the case for 
the third step (p=2) in the example of Fig. 20. The 3 
low-order bits of the second address generated for the 

25 first internal butterfly are 010. When the last bit (p- 
p-1— 1) of this number is taken away, the number 01 is 
obtained relating to coefficient W 1 . This case always 
corresponds to the last step of a transformation method 
comprising p-1 transformation steps. 

30 c) if p+Ku-2, the coefficient address corresponds 

to the number formed by the p+1 low-order bits of the 
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second address generated for this internal butterfly, 
followed by p-p-3 zero bits at the end of the number. 
This case is illustrated in Fig. 21. This figure 
represents the addresses relating to an internal 
5 butterfly for implementing the second transformation 
step (p=l) of a transformation method designed for 
processing a series of 32 real samples (N=32 and u=5) . 
This butterfly transforms the samples of rank 1, 2, 5, 
and 6 of the series of samples obtained at the end of 

10 the first step. The two low-order bits of the second 
generated address are 10 and if a zero (u-p-3=l) is 
added at the end of the number, the number 100 is 
obtained designating coefficient W 4 . 

Thus, both addresses produced by the address 

15 generator for one butterfly is used to address both 
four samples to be processed and the coefficient 
relating to the butterfly. 

Preferably, even address and odd address samples 
will be stored in two separate memories. Thus, it will 

20 be possible to read two input samples simultaneously, 
and it will be possible to write the resulting output 
samples simultaneously, which means saving process time 
for the series to be transformed. 

As mentioned before, for each calculation method 

25 operating according to a time interleaving algorithm 
there is a corresponding method operating according to 
a frequency interleaving algorithm. All that has to be 
done to obtain it is, on the one hand, to invert the 
series of transformation operations of the 

30 corresponding time interleaving method, and on the 
other hand, for each butterfly, to invert the 
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transformation operations themselves with respect to 
those of the corresponding method. 

In addition, the invention also relates to a 
method of calculating the fast Fourier transform or the 
5 inverse fast Fourier transform of a series of N complex 
samples X (n) conjugated by pairs, with N power of 2, 
operating according to a frequency interleaving 
algorithm. The series of N complex samples X(n) is 
represented by a series of N real samples y(n) defined 
10 as follows: 

y(0) - RetX(O) ] 

y(n) = Re [X ( (n+1 ) /2) ] for n being odd and 

different from N-l 
y(n) = Im[X(n/2)] for n being even and 

15 different from 0 

y (N-l) = Re [X (N/2) ] 

According to the invention, this method 
substantially comprises transformation steps for 
transforming input samples into output samples. Real 

20 samples y(n) are processed in a first transformation 
step and the last step provides a series of N output 
samples x(n) representative of the fast or the inverse 
fast Fourier transform of the sample series X(n). As 
for the time interleaving circuits, each transformation 

25 step is implemented by a set of butterflies with 
several inputs and several outputs . As the 
transformation is done according to an in place method, 
all the steps are performed by means of a single set of 
butterflies, the operating mode of which is modified in 

30 each transformation step. In each transformation step, 
input and output samples are stored in a storage 
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memory. When they have been transformed, the output 
samples of the same butterfly replace the corresponding 
input samples of the same rank in the storage memory. 

According to the invention, if the samples y(n) 
5 introduced in the first transformation step are 
classified in ascending order of index n, the output 
samples x (n) are provided in the last transformation 
step in bit-reversed order of index n. The output 
sequence of the x (n) samples can then be modified by a 
10 final step so as to classify them in ascending order of 
index n. 

An embodiment of such a transformation is 
represented in Fig. 22. It is inferred from the 
embodiment of Fig. 12 by inverting the functional 

15 arrangement of the transformation steps of Fig. 12 
(mirror image of what it is for time interleaving) . 
This embodiment allows the calculation of the inverse 
fast Fourier transform of a series of 16 real samples 
y(n) representative of a series of 16 complex samples 

20 X(n) conjugated by pairs. Coefficients W s are therefore 
of the type e 3<2,cs/N> . 

This embodiment comprises three transformation 
steps E p with 0<p<2 . In each transformation step E p , 
the butterflies are henceforth distributed among 2 P 

25 calculation blocks, these calculation blocks being 
sorted in each step according to an ascending rank q 
from 0 to 2 P -1. Each calculation block has a peripheral 
butterfly and N/2 p+2 -l internal butterflies. The y(n) 
samples are applied in the sequence of index n in the 

30 first step. 
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The ranks of the samples processed by the same 
butterfly are defined as follows: in step Ep, 
considering a peripheral butterfly belonging to the 
rank a calculation block, it transforms the input 
5 samples of rank 2 p - p a, 2 M_p a+2 lJ - p - 1 -l , 2^<x+2^~ 1 , 2 M - p a+2 p - 
p -l into output samples of the same rank, and, in step 
Ep, considering an internal rank x butterfly in the 
rank a calculation block, it transforms the input 
samples of rank 2 p ~ p oc+2x+l, 2 M-p a+2x+2, 2 VJ_p a+2 p_p -2x-3 , 2 M ~ 

10 p a+2 M_p -2x-2 into output samples of the same rank. 
Finally, the coefficient assigned to the rank x 
internal butterfly of the rank a calculation block in 
step Ep is equal to W 5 with 5 = (x+1) . 2 P . 

The real samples x (n) obtained at the end of the 

15 method are provided in bit-reversed order of index n. 

At the butterflies, the coefficients W s are of the 
type e J< 2,ts/N) an d inputs and outputs have been inverted 
with respect to the embodiment of Fig. 12. 
Consequently, the operations performed by the 

20 butterflies of this embodiment are different from those 
performed by the butterflies of Fig. 12. A butterfly 
design associated with the embodiment of Fig. 22 is 
represented in Fig. 23. 

Just like the design associated with Fig. 12, it 

25 has four data inputs and four data outputs as well as 
four additional inputs, respectively primary mode MP, 
secondary mode MS, permutation PERM, and coefficient 
COEF inputs. The primary mode signal is 0 for an 
peripheral butterfly and 1 for an internal butterfly. 

30 The permutation signal is 0 for even values of rank q 
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and 1 for odd values. Finally, the secondary mode 
signal is 1 if the peripheral butterfly is used for 
implementing the first step, and otherwise 0. 

The calculations performed by this butterfly are 
5 the following ones (W s =A+j.B) : 

1) if the primary mode and secondary mode signals 
are 0 : 

si = (el + e2) /2 
s2 = (el - e2) /2 
10 s3 = (e4 - e3) /2 

s4 = (e3 + e4) /2 

2) if the primary mode signal is 0 and the 
secondary mode signal is 1 : 

si = [ (el+e4) /2 + e2] /2 
15 s2 = [ (el + e4) /2-e2] /2 

s3 = - [e3- (el-e4) /2] /2 
s4 = [e3+ (el-e4) /2] /2 

3) if the primary mode signal is 1 and the 
permutation signal is 0: 

20 si = (el+e3)/2 

s2 = (e2-e4)/2 

s3 = [ (el-e3) .A - (e2+e4).B]/2 
s4 = [-(el-e3).B + (e2+e4).A]/2 

4) if the primary mode signal is 1 and the 
25 permutation signal is 1 : 

si = [ (el-e3) .A - (e2+e4).B]/2 
s2 = [-(el-e3).B + (e2+e4).A]/2 
s3 = (el+e3) /2 
s4 = (e2-e4)/2 

30 This design is inferred from the design of Fig. 13 

by inverting values el, e2, e3, e4 and values si, s2, 
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s3, s4, and by replacing B with -B because coefficient 
W s is now of the type e j(2ns/N) . As A and B represent the 
cosine and sine of the same number, they have become 
A 2 +B 2 =l . The expressions of si, s2, s3, s4 are thus 
5 simplified. 

According to the invention, addresses are 
furthermore generated for sample addressing. Provision 
is made for generating two binary addresses of u bits 
per butterfly, each binary address corresponding to the 

10 rank of a butterfly input sample. The addresses of the 
other two samples to be applied to the butterfly inputs 
are obtained by inverting the u-p low-order bits of the 
first two addresses. 

In the same way as for the transformation methods 

15 operating according to time interleaving, the two 
binary addresses produced are consecutive for an 
internal butterfly. For a peripheral butterfly, the u-p 
low-order bits of the first generated address are equal 
to 0, and the u-p low-order bits of the second address 

20 form a number equal to N/2 P+1 -1 . By way of example, the 
addresses produced for the transformation circuit of 
Fig. 22 are gathered in Fig. 24. 

Advantageously, provision can be made for storing 
even address samples and odd address samples in two 

25 separate memories in order to reduce processing time of 
the transformation operation. 

Finally, the addresses generated for addressing 
the samples are also used for addressing coefficients 
W s . The value of the parameter s is used for addressing 

30 the corresponding coefficient W s . In this embodiment, 
the parameter s is equal to: 
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if u-p-l=u-2, the number formed by the u-p-1 
low-order bits of the second address produced for said 
internal butterfly, 

if u-p-Ku-2, the number formed by the u-p-1 
5 low-order bits of the second address produced for said 
internal butterfly, followed by p-1 zero bits at the 
end of the number, 

if u-p-l>u-2, the number formed by the u-p-1 
low-order bits of the second address produced for said 
10 internal butterfly, minus its p+1 low-order bits. This 
case corresponds to the first step (p=0) of the 
transformation methods operating according to a 
frequency interleaving algorithm. 
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CLAIMS 



1. A method of calculating the fast Fourier 
transform or the inverse fast Fourier transform of a 
5 digital signal defined by a series of N real starting 
samples x (n) , with N a power of two and n e [0..N-1], 
comprising successive transformation steps (2) for 
transforming input samples into output samples, all the 
transformation steps being performed by means of a 

10 single set of butterflies with several inputs and 
several outputs, the operating mode of which is 
modified selectively in each transformation step, the 
input and output samples of each transformation step 
being stored in a storage memory, a series of N output 

15 samples y(n) representative of the fast Fourier 
transform or the inverse fast Fourier transform of the 
output samples x (n) being provided in the last 
transformation step, 

characterized in that output samples y(n) are 

20 real, 

and in that the output samples of a butterfly 
replace the corresponding input samples of the same 
rank in the storage memory, so that, if the starting 
samples x (n) processed in the first transformation step 
25 are classified in bit-reversed order of their index n, 
output samples y(n) are provided in the last 
transformation step in ascending order of index n, 
these output samples being defined by the following 
relations : 
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y (0) = Re [X (0) ] 

y(n) = Re [X ( (n+1) /2) ] for n being odd and 
different from N-l 

y(n) = Im[X(n/2)] for n being even and 

different from 0 

y (N-l) = Re [X (N/2) ] 
where samples X (n) , with n e [0..N-1], designate the 
complex samples of the series corresponding to the fast 
or inverse fast Fourier transform of the starting 
sample series x(n) . 

2. A method of calculating the fast Fourier 
transform or the inverse fast Fourier transform of a 
digital signal defined by a series of N complex samples 
X(n) conjugated by pairs represented by a series of N 
real starting samples y(n), with N power of two and n 
e [0.. N-l], the starting samples y(n) being defined as 
follows: 

y (0) = Re [X (0) ] 

y(n) = Re [X ( (n+1) /2) ] for n being odd and 
different from N-l 

y(n) = Im[X(n/2)] for n being even and 

different from 0 

y (N-l) - Re [X (N/2) ] 
this calculation method comprising successive 
transformation steps for transforming input samples 
into output samples, a series of N output samples x(n) 
representative of this fast or inverse fast Fourier 
transform being provided in the last transformation 
step, all the transformation steps being performed by 
means of a single set of butterflies with several 
inputs and several outputs, the operating mode of which 
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is modified selectively in each transformation step, 
the input and output samples of each transformation 
step being stored in a storage memory, 

characterized in that output samples x (n) are 

real, 

and in that the output samples of a butterfly 
replace the corresponding input samples of the same 
rank in the storage memory, so that, if the starting 
samples y(n) processed in the first transformation step 
are classified in ascending order of index n, the 
output samples x (n) are provided in the last 
transformation step in bit-reversed order of index n. 

3. The calculation method according to claim 1 or 
2, characterized in that, in each transformation step, 
each butterfly transforms input sample pairs, the ranks 
of the input samples of the same pair within the series 
of input samples of said transformation step being 
symmetrical with respect to a center between the end 
rank values of the input samples transformed by said 
butterfly. 

4. The calculation method according to claim 3, 
characterized in that it comprises u-1 transformation 
steps E p with u=log 2 (N) and p e [0..U-2]. 

5. The calculation method according to claim 4, in 
turn dependent on claim 3, in turn dependent on claim 
1, characterized in further comprising: 

a preliminary step of modifying the sequence of 
the starting samples x(n) ranked in ascending order of 
index n and showing them in bit-reversed order of index 
n in the first transformation step, and 



- a final step of processing the series of output 
samples y(n) and providing a series of N complex 
conjugated samples X (n) corresponding to the fast or 
the inverse fast Fourier transform of the series of 
starting samples x(n) . 

6. The calculation method of claim 4, in turn 
dependent on claim 3, in turn dependent on claim 1, or 
according to claim 5, characterized in that, in each 
transformation step E p/ butterflies are distributed 
among N/2 P+2 calculation blocks, 

in that each calculation block has a peripheral 
butterfly and/or 2 P -1 internal butterflies, 

in that the peripheral butterfly of the rank a 
calculation block in transformation step Ep transforms 
the input samples of rank 2 p+2 a, 2 p+2 a+2 p+1 - 1 , 2 p+2 a+2 p+1 , 
2 p+2 a+2 p+2 -l into output samples of the same rank, 

and in that the internal rank x butterfly of the 
rank a calculation block in transformation step Ep 
transforms the input samples of rank 2 p+2 a+2x+l, 
2 p+2 a+2x+2, 2 p+2 a+2 p+2 -2x-3, 2 p+2 a+2 p+2 -2x-2 into output 
samples of the same rank, with P>1 . 

7. The calculation method according to claim 6, 
characterized in that each butterfly is assigned a 
coefficient W s , whereon the calculation inside the 
butterfly is based, said coefficient being equal to e~ 
j(2n s /N) with s G [0.. N/4-1] for a fast Fourier transform 
and is equal to e j(27ts/N) with s <= [0.. N/4-1] for an 
inverse fast Fourier transform. 

8. Calculation method according to claim 7, 
characterized in that the internal rank x butterfly of 



40 



the rank a calculation block in transformation step Ep 
is assigned coefficient W 8 with 5 = (x+1) . (N/2 P+2 ) . 

9. The calculation method according to claim 8, 
characterized in that the butterflies for implementing 

5 the transformation steps are all of the same type and 
have 

- four inputs for receiving input samples and 
four outputs for providing output samples, 

- four additional inputs, respectively primary 
10 mode, secondary mode, permutation, and coefficient 

inputs, 

in order to selectively apply different 
transformation operations to the input samples, each 
operation being determined by the values assigned to 
15 the primary mode, secondary mode, permutation signals, 
and a coefficient admitted on said corresponding 
additional inputs. 

10. The calculation method according to claim 9, 
characterized in that, for each butterfly, the primary 

20 mode signal is 0 for a peripheral butterfly and 1 for 
an internal butterfly, 

in that the permutation signal is 0 for the even 
rank calculation blocks, including rank 0, and 1 for 
the other ones. 

25 11. The calculation method according to claim 10, 

characterized in that, in transformation step Ep, each 
calculation block comprises one peripheral butterfly 
and 2 P -1 internal butterflies. 

12. The calculation method according to claim 11, 

30 characterized in that the secondary mode signal is 1 if 
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the peripheral butterfly is used for the last 
transformation step, and otherwise 0. 

13. The calculation method according to claim 12, 
characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W S =A+ j . B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if the primary mode and secondary mode signals 
are 0: si = el + e2 

s2 = el - e2 
s3 = e4 - e3 
s4 = e3 + e4 

2) if the primary mode signal is 0 and the 
secondary mode signal is 1: 

si - el + e2 + e3 + e4 
s2 = el - e2 
s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

3) if the primary mode signal is 1 and the 
permutation signal is 0: 

si = el + A.e3 - B.e4 
s2 = e2 + B.e3 + A.e4 
s3 = el - A.e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 

4) if the primary mode signal is 1 and the 
permutation signal is 1: 

si = el - A.e3 + B.e4 
s2 = -e2 + B.e3 + A.e4 
s3 = el - A.e3 - B.e4 
s4 = e2 + B.e3 + A.e4 
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14. The calculation method according to claim 10, 
characterized in that, in transformation step E p , each 
calculation block comprises: 

2 P -1 internal butterflies and a peripheral 
5 butterfly for the even values of index p as well as for 
the last transformation step if u is even, and 
2 P -1 internal butterflies, otherwise. 

15. The calculation method according to claim 13, 
characterized in that the secondary mode signal is 1 if 

10 the peripheral butterfly is used for the last 
transformation step with u being odd, and otherwise 0. 

16. The calculation method according to claim 15, 
characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W s =A+j.B, the 

15 butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if primary mode, secondary mode and permutation 
signals are 0: 

si = el + e2 + e3 + e4 
20 s2 = el - e2 

s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

2) if the primary mode signal is 0 and the 
secondary mode signal is 1 : 

25 si = el + e4 

s2 = e2 
s3 = e3 
s4 = el - e4 

3) if the primary mode signal is 0 and the 
30 permutation signal is 1 : 



4) if the primary mode signal is 1 and the 
permutation signal is 0: 

si = el + A.e3 - B.e4 
s2 - e2 + B.e3 + A.e4 
s3 = el - A.e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 

5) if the primary mode signal is 1 and the 
permutation signal is 1: 

si = el - A.e3 + B.e4 
s2 = -e2 + B.e3 + A.e4 
s3 = el + A.e3 - B.e4 
s4 = e2 + B.e3 + A.e4 

17. The calculation method according to claim 10, 
characterized in that, in transformation step E p , each 
calculation block comprises: 

2 P -1 internal butterflies and a peripheral 
butterfly for the even values of index p, and 
2 P -1 internal butterflies, otherwise. 

18. The calculation method according to claim 17, 
characterized in that the secondary mode signal is 1 if 
the peripheral butterfly is used for the first 
transformation step with u being even, and otherwise 0. 

19. The calculation method according to claim 18, 
characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W s =A+j.B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 
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1) if the primary mode signal is 0 and the 
secondary mode signal is 1: 

si = el + e2 
s2 = el - e2 
s3 - e4 - e3 
s4 = e3 + e4 

2) if primary mode, secondary mode and permutation 
signals are 0: 

si = el + e2 + e3 + e4 
s2 = el - e2 
s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

3) if the primary mode and secondary mode signals 
are 0 and the permutation signal is 1 : 

si = (e3 + e4) - (el + e2) 

s2 = el - e2 

s3 = e4 - e3 

s4 = el + e2 + e3 + e4 

4) if the primary mode signal is 1 and the 
permutation signal is 0: 

si = el + A.e3 - B.e4 
s2 = e2 + B.e3 + A.e4 
s3 = el - A.e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 

5) if the primary mode signal is 1 and the 
permutation signal is 1 : 

si = el - A.e3 + B.e4 
s2 = -e2 + B.e3 + A.e4 
s3 = el + A.e3 - B.e4 
s4 = e2 + B.e3 + A.e4 
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20. The calculation method according to claim 8, 
characterized in that the butterflies for implementing 
the transformation steps are all of the same type and 
have 

5 - four inputs for receiving input samples and 

four outputs for providing output samples, 

four additional inputs, respectively primary 
mode, secondary mode, permutation, and coefficient 
inputs, 

10 in order to selectively apply different 

transformation operations to the input samples, each 
operation being determined by the values assigned to 
the primary mode, secondary mode, permutation signals, 
and a coefficient admitted on said corresponding 

15 additional inputs, 

and in that the final step furthermore performs an 
addition and subtraction between the first and the last 
output sample provided in the last transformation step. 

21. The calculation method according to claim 20, 
20 characterized in that, in transformation step E p , each 

calculation block comprises one peripheral butterfly 
and 2 P -1 internal butterflies. 

22. The calculation method according to claim 21, 
characterized in that, for four input samples el, e2, 

25 e3, and e4, and for a complex coefficient W s =A+j.B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if the primary mode signal is 0: 
si = el + e2 
30 s2 = el - e2 
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2) if the primary mode signal is 1 and the 
permutation signal is 0: 

si = el + A.e3 - B.e4 
5 s2 = e2 + B.e3 + A.e4 

s3 = el - A. e3 + B . e4 
s4 = -e2 + B.e3 + A.e4 

3) if the primary signal is 1 and the permutation 
signal is 1: 

10 si = el - A.e3 + B.e4 

s2 = -e2 + B.e3 + A.e4 
s3 = el + A.e3 -B.e4 
s4 = e2 + B.e3 + A.e4 
23. The calculation method according to claim 9 or 
15 20, characterized in that the first and second binary 
addresses of u bits are generated for each butterfly, 
each binary address corresponding to the rank of an 
input sample of said butterfly and the second binary 
address being greater than the first binary address. 
20 24. The calculation method according to claim 23, 

characterized in that said first and second binary 
addresses are consecutive and an internal butterfly is 
involved . 

25. The calculation method according to claim 23 
25 or 24, characterized in that, if a peripheral butterfly 

is involved, the p+2 low-order bits of the first 
address are equal to 0, and the p+2 low-order bits of 
the second address form a number equal to 2 P+1 -1. 

26. The calculation method according to claim 24 
30 or 25, characterized in that the address of the two 

other samples to be applied to the inputs of the 
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butterfly, be they peripheral or internal, are obtained 
by inverting the (p+2) low-order bits of said first and 
second produced addresses. 

27. The calculation method according to claim 26, 
characterized in that even-numbered address samples and 
odd-numbered address samples are stored in two separate 
memories . 

28. The calculation method according to claim 25, 
characterized in that the value of the parameter s of 
the coefficient W s assigned to an internal butterfly in 
transformation step E p is coded by u-2 bits, and is: 

- if p+l=u-2, the number formed by the p+1 
low-order bits of the second binary address produced 
for said internal butterfly, 

- if p+Ku-2, the number formed by the p+1 
low-order bits of the second binary address produced 
for said internal butterfly, followed by u-p-3 zero 
bits at the end of the number, 

if p+l>p-2, the number formed by the p+1 low- 
order bits of the second binary address produced for 
said internal butterfly, minus its p-p-1 low-order 
bits. 

29. The calculation method according to claim 4, 
in turn dependent on claim 3, in turn dependent on 
claim 2, characterized in that in each transformation 
step E p , the butterflies are distributed among 2 P 
calculation blocks, 

in that each calculation block comprises one 
peripheral butterfly and N/2 p+2 -l internal butterflies, 

in that the peripheral butterfly of the rank a 
calculation block in transformation step Ep transforms 
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the input samples of rank 2 VJ_p a, 2 p_p a+ 2 M-P_1 - 1 , 2 p_p ot+ 2 VJ_P " 
l , 2 M-P a+2 1J ~ P -1 into output samples of the same rank, 

and in that the internal rank x butterfly of the 
rank a calculation block in transformation step Ep 
5 transforms the input samples of rank 2 M ~ p a+2x+l, 2 U ~ 
p a+2x+2, 2 M-p a+2 p ~ p " 1 -2x-3, 2 M ~ p a+2 M_p -2x-2 into output 
samples of the same rank. 

30. The calculation method according to claim 29, 
characterized in further comprising a final step of 

10 modifying the sequence of the output samples provided 
in the last transformation step and classifying them in 
ascending order of index n. 

31. The calculation method according to claim 29 
or 30, characterized in that each butterfly is assigned 

15 a coefficient W s , whereon the calculation inside the 
butterfly is based, said coefficient being equal to e" 
j(2ns/N) with s e [0.. N/4-1] for a fast Fourier transform 
and is equal to e j(2ns/N) with s e [0.. N/4-1] for an 
inverse fast Fourier transform. 

20 32. Calculation method according to claim 31, 

characterized in that the internal rank x butterfly of 
the rank a calculation block in transformation step Ep 
is assigned coefficient W s with 5 = (x+l).2 p . 

33. The calculation method according to claim 32, 

25 characterized in that the butterflies for implementing 
the transformation steps are all of the same type and 
have 

- four inputs for receiving input samples and 
four outputs for providing output samples, 
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- four additional inputs, respectively primary 
mode, secondary mode, permutation, and coefficient 
inputs, 

in order to selectively apply different 
5 transformation operations to the input samples, each 
operation being determined by the values assigned to 
the primary mode, secondary mode, permutation signals, 
and a coefficient admitted on said corresponding 
additional inputs. 
10 34. The calculation method according to claim 33, 

characterized in that, for each butterfly, the primary 
mode signal is 0 for a peripheral butterfly and 1 for 
an internal butterfly, 

in that the permutation signal is 0 for the even 
15 rank calculation blocks, including rank 0, and 1 for 
the odd values. 

35. The calculation method according to claim 31 
or 34, characterized in that the secondary mode signal 
is 1 if the butterfly, be it peripheral or internal, is 

20 used for the first transformation step, and otherwise 
0. 

36. The calculation method according to claim 35, 
characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W s =A+j.B, the 

25 butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if the primary mode and secondary mode signals 
are 0: 

si = (el + e2)/2 
30 s2 = (el - e2)/2 

s3 = (e4 - e3)/2 
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s4 = (e3 + e4) /2 

2) if the primary mode signal is 0 and the 
secondary mode signal is 1: 

si = [ (el + e4) /2-e2] /2 
s2 = [ (el+e4) /2-e2] /2 
s3 = [e3- (el-e4) /2] /2 
s4 - [e3+ (el + e4) /2] /2 

3) if the primary mode signal is 1 and the 
permutation signal is 0: 

si = (el+e3)/2 
s2 = (e2+e4)/2 

s3 = [ (el-e3) .A - (e2 + e4) .B]/2 
s4 = [-(el-e3).B + (e2+e4).A]/2 

4) if the primary mode signal is 1 and the 
permutation signal is 1: 

si = [(el-e3).A- (e2+e4).B]/2 
s2 = [-(el-e3).B + (e2+e4).A]/2 
s3 = (el+e3)/2 
s4 = (e2-e4)/2 

37. The calculation method according to claim 33, 
characterized in that the first and second binary 
addresses of u bits are generated for each butterfly, 
each binary address corresponding to the rank of an 
input sample of said butterfly and the second binary 
address being greater than the first binary address. 

38. The calculation method according to claim 37, 
characterized in that said first and second binary 
addresses are consecutive and an internal butterfly is 
involved. 

39. The calculation method according to claim 37 
or 38, characterized in that, if a peripheral butterfly 
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is involved, the u-p low-order bits of the first 
address are equal to 0, and the u-p low-order bits of 
the second address form a number equal to N/2 P+1 -1. 

40. The calculation method according to claim 38 
5 or 39, characterized in that the address of the two 

other samples to be applied to the inputs of the 
butterfly are obtained by inverting the u-p low-order 
bits of both produced addresses. 

41. The calculation method according to claim 40, 
10 characterized in that even-numbered address samples and 

odd-numbered address samples are stored in two separate 
memories. 

42. The calculation method according to claim 41, 
characterized in that the value of the parameter s of 

15 the coefficient W s assigned to an internal butterfly in 

transformation step E p is coded by u-2 bits, and is: 

- if p-p-l=u-2 , the number formed by the u-p- 1 

low-order bits of the second address produced for said 

internal butterfly, 
20 - if u-p-Ku-2, the number formed by the u-p-1 

low-order bits of the second address produced for said 

internal butterfly, followed by p-1 zero bits at the 

end of the number, 

if u-p-l>u-2, the number formed by the u-p-1 
25 low-order bits of the second address produced for said 

internal butterfly, minus its p+1 low-order bits. 



ABSTRACT 

METHOD OF CALCULATING THE FAST FOURIER TRANSFORM AND 
THE INVERSE FAST FOURIER TRANSFORM 

5 

The invention relates to a method of calculating 
the fast Fourier transform or the inverse fast Fourier 
transform of a series of N real samples x (n) , with N 
power of two, operating according to a time 

10 interleaving algorithm and providing the sample series 
X (n) in ascending order of index n and using limited 
calculation and storage means. The invention also 
relates to a method of calculating the fast Fourier 
transform or the inverse fast Fourier transform of a 

15 series of N conjugated complex samples X (n) , with N 
power of two, operating according to a frequency 
interleaving algorithm. 
Fig. 8 

Application: image or acoustic signal processing, 
20 multicarrier modulation. 
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[x(0) ,x(l) ,x(2) ,x(3) ,x(4) ,x(5) ,x(6) ,x(7) , 
x(8) ,x(9) ,x(10) ,x(ll) ,x(12) ,x(13) ,x(14) f x(15) ] 



PRELIMINARY STEP 



[x(0) ( x(8),x(4),x(12),x(2),x(10),x(6),x(14) , 
x(l),x(9),x(5),x(13),x(3),x(ll),x(7),x(15)] 



TRANSFORMATION STEPS 



CY(0) ,y(i) ,y(2) ,y(3) ,y(4) ,y(5) ,y(6) ,y(7) , 
Y(8) ,y(9) ,y(10) ,y(ll) ,y(12) ,y(13) ,y(14) ,y(15) ] 



FINAL STEP 



[X(0) ,X(1) ,X(2) ,X(3) ,X(4) ,X(5) ,X(6) ,X(7) , 
X(8) ,X(9) ,X(10) ,X(11) ,X(12) ,X(13) ,X(14) ,X(15) ] 



FIG. 6 
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|x(2)+x(10)=A^4T 



|x(6)+x(14)=>W6) 



|B B (13) 



I -S,(13) 
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7> 

A«(4)-t-A«(6)=B R (4) 



. A«(4)-Aa(6)=B R (6) 



FIG.7A 



B R (9)+3.B R (1 3)+3~B?(1 3)=C r (9)~ 
B,(9)-3.B r (1 3)+3.B,(1 3)=C,(9) 



, B R (9)-3.B„(13)-3.B,(13)=C a (11) 
-B,(9)-3.B W (1 3)+3.B,(1 3)=C,(1 1 ) 



- I -[A r (4)-Ak(6)]=-B r (6) 



lA^f4)-FA«(6)=B R (4) 




8 R (9)-3.B,(13)-3.Bi(13)=C B (11) 
-B,(9)-3.B R (13)»3.B,(13)=C,(1 1) 



■B R (9)+3.B R (13)+3.B|(13)=C R (9) | 



B,(9)-3.B R (1 3)+3.B,(1 3)=C,(9) 



FIG.7B 




FIG. 9 
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I haroby dalm ffvo benefit undar 35 U.S.C. HB(o) of any United States provisional application^) llsied twlow, 

— , Am i „ au „ ■ rr. — TTT — r: = — D AddlUonat provisional application 

(AeplicftHon Number) (MoniWDay/Year Rlad) numben are lltlad on a 

supplemental priority data sha^l 



(Application Numbar) (ManlfwOay/Yoar filed) 



PTO/SB/02B ^llaehed hereto. 



OCT 17 2002 09=30 



NILLES NILLES 



ii'ti a, it ■!■.'„ ;B i! ; 41^ , 2 7fe ''ifefe '"' ! ' !i r*; ! 03 

1 1 !> & 4 



DECLARATION - Utility or Design Patent Application 

Inventors: AJj JALALI, Pierre LeRAY. Dominique LACROIX 
Serial No. 09/581.272 

I hereby claim the benefit under 35 U.5.C 120 of any United Stales appllaWonjg). or 365(c) of any PCT Intemullanal application 
designallnc; the UnliM States, of America, listed below and. Insofar as lha subject matter or each of the claims of w» application is 
not diadosad In (he prior Untied Slates or PCT International application In tho manner pnwJd&d by the first paragraph or 35 U.S.C. 
112, I acknowledge the truly to disclose Information *hieb Is material lo the patentability ai defined, in 37 CFR 1.56 wnich boearnq 
available between Via tiling data of the prior application and tna national or PCT International filing data of this application. 
U.S. Parent Application or PCT Parent Applications) 



PCT/FR9tV0283a 



■ 2000 



(MontMDay/Year Filed) 



(Patent Number (If applicable)) 



(NumDer) (Montn/Day/toor Filed) (Patent Number (if applicable)) 

).S. or PCT taxational npptaUort numbers are lis tad on a supplemental priority data sheet PTO/SBr32S attached 



James 6. Nll(a«, Reg. No. 16,603 
Andraw J- Nilloa, Rag. No, 2AJ»&~^ 
Jay O. Dursl. Reg. No^4i723 
Usa M. Getirko, Reg. TtoTSRHBtTl 
Tftaddau* C. EUnkowsKI, Reg. No. 45,522 

Direct all talaphone calls to Jamoa E, Mllaa at ajtaphpna nur 
Direct all corrospondoncc lo: Jamas E. Nlllas 

TRIlDESTWCEIs. S.C. 
"eKSrSKl 

"777 Eoat Wbcortali 

Uln 53202-53*5 



Matthew C. Loppnow. Row- No. «,314 
Stopnon Mlcnaoi Patton. Rag. >Jo. 3fl,23S 
jBromo D. Dreplak, Reg. No. 31.011 
Lisa A. BrzycW. Rafl. No. 40.926 



(41*) 27M977. faeotmilo number (*14) 276-09H2. 



I horooy dodara that all statements mado herein of my own kncMrtodgo lira true ana that all statement* mado <n IntbrmaUgn and 
boliar an believed to ba true: and further thai these stalomonu vera made with the knowledge that vwflfrtif (also statements and tea 
OKa ao made afo punishable Dj fine or Imprisonment or both, under 18 U.S.C. 1001 and thai such willful falsa statements may 
Jeopardize the validity of ma application or any palanl Issued thereon. 

Full nomo of Sole or First Inventor: 0 a petition nas been filed to this unsfanad Inventor 

Gtvan Name (first & middle fjf any]) 4 Family NarrWSvmame: All JALAU 



Inwonjor'a Signature 



ilOence (city, state, country): Rennea. Franca 
Post Offlca Address: 32 rue Mlrabeciu 

Bailment K No. 234 
(dry, «b(a. Op. country): 35700 Rennes. Fram 



OCT-17-2002 09=30 



N1LLES NILLES 



DECLARATION - Utility or Design Patent Application 

Inventors: AH JALALI. Pierre LeRAY. Dominique LACROIX 
Serial No. 09/581.272 

Full nama of Second ln«an«or, if any: D A pwlidon Ubb boon filed for Uiis unsigned Invonlor 

GtVon Nnmn (finrt & middle (If wyl) & Family Name/Surname: Pierre LERAY 

Invent. SX,natu«: j^m^J^SSAE-.. ^f^^^ t t . Date */£oOO 

/T>» Residence (dry, efcale. country): Uffro, Franca f~M^ ^J^^ ' CltlzeftsNp; French 

(/ Post Offlai Address: LaGrass* Roche 

(dV. *»te. zip. cauniry): SSKO Uffoe. Fninaa 

Full nama of Thlfd Invantor, If any: Q A petition l\a* Been filed for <Ma unsigned Invanlor 

Given Name (Aral & middle [if any]) & Family Name/Sum amn: Domlnlauft LACROIX 



5' 



inventor-, Sfaasmro; £&m&dAM± £Jit#Q 'l X ( =F%3 && $£^ 1A 1 1 I *20OO 

Residence [dry, etata. cpuniry) ^ Ratines, F rance /"^X,— — ' > Cllizenshlp: French 

0^ Post Oflien AddrQ»: w > Cwu i u Ala i n r ci w it 3^ S^awc l^xift 6p>£x^«V 



Page 3 of 3 



